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In this paper, a nonlinear trajectory control algorithm of rendezvous with maneuvering 
target spacecraft is presented. The disturbance forces on the chaser and target spacecraft 
and the thrust forces on the chaser spacecraft are considered in the analysis. The control 
algorithm developed in this paper uses the relative distance and relative velocity between the 
target and chaser spacecraft as the inputs. A general formula of reference relative trajectory 
of the chaser spacecraft to the target spacecraft is developed and applied to four different 
proximity maneuvers, which are in-track circling, cross-track circling, in-track spiral 
rendezvous and cross-track spiral rendezvous. The closed-loop differential equations of the 
proximity relative motion with the control algorithm are derived. It is proven in the paper 
that the tracking errors between the commanded relative trajectory and the actual relative 
trajectory are bounded within a constant region determined by the control gains. The 
prediction of the tracking errors is obtained. Design examples are provided to show the 
implementation of the control algorithm. The simulation results show that the actual relative 
trajectory tracks the commanded relative trajectory tightly. The predicted tracking errors 
match those calculated in the simulation results. The control algorithm developed in this 
paper can also be applied to interception of maneuver target spacecraft and relative 
trajectory control of spacecraft formation flying. 


I. Introduction 

T RAJECTORY control of spacecraft rendezvous has been studied for many years. For the most of the studies, the 
target spacecraft doesn’t maneuver by firing thrusters during the rendezvous [1-5]. Using the linearized 
Tschauner-FIempel equations, spacecraft rendezvous with target spacecraft in an arbitrary elliptical orbit was 
analyzed [1]. The constrained regulation problem is solved based on a parametric Lyapunov differential equation 
approach. The disturbance forces on the chaser and target spacecraft and the thrust forces on the chaser spacecraft 
are not considered in the analysis in [1], The six-degree-of-freedom elliptical orbit rendezvous and docking problem 
were investigated using an adaptive discrete-time robust nonlinear filtering technique [2]. However, the disturbance 
and thrust forces on the target spacecraft are neglected in the relative motion dynamics in [2]. 

In this paper, a nonlinear trajectory control algorithm of rendezvous with maneuvering target spacecraft is 
presented. The relative motion equation of the maneuvering target spacecraft and chaser spacecraft is developed. 
The acceleration due to thrust forces on the target spacecraft for maneuvers is included in the equation. The 
disturbance forces, such as aerodynamic drag and solar pressure are also considered in the equation. The relative 
motion is governed by a nonlinear time variant differential equation. 

The target spacecraft orbit variation, due to its maneuvers, imposes a challenge for the control system design. 
The control algorithm developed in this paper uses the relative distance and relative velocity between the target and 
chaser spacecraft as the inputs. The relative distance and relative velocity can be measured by a lidar sensor. A 
reference relative trajectory of the chaser spacecraft to the target spacecraft are formed based on the required 
rendezvous, such as approaching the target spacecraft and circling it or approaching the target spirally and docking 
with it eventually. The reference relative trajectory can be in the in-track plan of the target spacecraft or the cross- 
track plan. The thrust command output from the control system is calculated based on the differences of the relative 
distance and the reference relative trajectory as well as the measured relative velocity. 

Several proximity and rendezvous maneuvers are realized using the control algorithm described in this paper. 
The first proximity maneuver is in-track circling, in which the chaser spacecraft circulates the target spacecraft at a 
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circular orbit relative to the target with a constant rate. The relative circular orbit is in the target spacecraft’s in- 
track plan. Similarly the cross-track circling is to circulate the target spacecraft in the target spacecraft’s cross-track 
plan. In the next two rendezvous maneuvers the chaser spacecraft approaches the target spacecraft along a spiral 
trajectory in the target spacecraft’s in-track or cross-track plans. Design examples and simulation results are 
presented at the end of the paper to demonstrate the performance of the control algorithm. The control algorithm 
developed in this paper can be applied to interception of maneuver target spacecraft and relative trajectory control of 
spacecraft formation flying also. 


II. Proximity Relative Motion Model 

The proximity relative motion model of chaser spacecraft relative to target spacecraft has to be established 
before the rendezvous trajectory control algorithms can be developed. In Figure 1, the target spacecraft orbit frame 

X ot ■ y m > Z oi an d ^e chaser spacecraft orbit frame X oc ,y oc , Z oc are depicted. The trajectory dynamics for the 
target spacecraft in the inertial frame is 


R 

R - - /t— j +a dt + a , 


(i) 


R is the distance vector from the center of gravity to the target spacecraft and jU is the gravitational parameter of 


the planet. a dt is the target spacecraft acceleration with respect to (wrt) the inertial frame due to the disturbance 


forces, such as aerodynamic drag, solar pressure. a t is the target spacecraft acceleration wrt the inertial frame due to 
the spacecraft thrust forces. 



Figure 1. Spacecraft redezvous geometry and coordinates. 


The equation of motion of the chaser spacecraft in the inertial frame is given by 

R + r 


R, = R + r = —/u - 


i? + r 


3 + a dc + Q c 


( 2 ) 


r is the relative distance vector from the target spacecraft to the chaser spacecraft. a dc is the target spacecraft 


acceleration wrt the inertial frame due to the disturbance forces. O c is the target spacecraft acceleration wrt the 

inertial frame due to the thrust forces. Using Eq.(l) and Eq.(2), the differential equation of the relative motion can be 
obtained as 
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r=M 

Note that the above equation has simple form because the differential equation is written in the inertial frame X J , 

Yj , Zj . In general, the disturbance forces and the thrust forces are calculated in the spacecraft body frame. The 
accelerations in the inertial frame can be obtained by transform the disturbance forces and the thrust forces from the 
body frame to the inertial frame as 


R 

Lk | 3 


R + r 
I R + H 3 


(3) 


A L is 


a dt - 

_ 4, F * 

(4) 

m, 



A!T 

(5) 

a, = 

. t 


m, 


a dc 


(6) 

m c 


a, = 

_ A X 

(7) 


m c 


target 

spacecraft body frame to the inertial frame. A bc 

is the 


transformation matrix from the chaser spacecraft body frame to the inertial frame. F dt is the disturbance forces on 
the target spacecraft wrt the target spacecraft body frame. m t and m are the target and chaser spacecraft mass 
respectively. T j is the target spacecraft thrust forces wrt its own body frame. However, T c is the chaser spacecraft 


thrust forces wrt the target spacecraft orbit frame X ot , y t , Z ot because the rendezvous trajectory control algorithms 
are formed in the target spacecraft orbit frame as described in the next section. The chaser spacecraft thrust forces 
wrt its own body frame T bc can be calculated by 

T bc = A b 0 X ( 8 ) 

Eq.(3) describes the relative motion between the target and the chaser spacecraft. No simplification and 
approximation are made in Eq.(3). The equation includes all the forces acting on the target and the chaser spacecraft. 
The target spacecraft can perform any maneuver by firing the thrusters. 


Ill, Trajectory Control Algorithms and Performance 

This section develops the trajectory control algorithm to control the chaser spacecraft relative motion to the target 
spacecraft. The control algorithm calculates the chaser spacecraft thrust forces wrt the target spacecraft orbit frame 
T„ as 


T c = m c \ K r {r cmd - A i' r )+ K v {r cmd - A °‘r)+ r cmd } (9) 

K r and K v are control gains, which can be matrix or scalar, r is the relative distance vector from the target 


spacecraft to the chaser spacecraft expressed in the inertial frame. A ° is the transformation matrix from the inertial 

frame to the target spacecraft orbit frame, which transforms r to the orbit frame, r is the relative velocity, which is 
the derivative of the relative distance vector wrt the inertial frame. The relative distance and relative velocity can be 
measured by a lidar sensor. r cmd is the commanded reference relative distance, which is formed based on the type of 
the relative motions, such as approaching the target spacecraft and circling it or approaching the target spirally and 
docking with it eventually. r cmd is the commanded relative velocity, which is the derivative of the commanded 


relative distance vector wrt the target orbit frame. r cmd is the derivative of r cmd wrt the target orbit frame. 

r cmd can be determined in the target orbit frame by the azimuth angle (X and the elevation angle (5 as shown in 
Figure 2. The commanded reference relative distance is given by 
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P.cn.AP" 


r ycmd (0 

= 

/w( 0, 

V 


' r m (t) cos a (0 cos /3(tp 

/(0 = r ycmd {t ) = r m (t) sin a(t) cos Pit) (10) 

f m (t) is the magnitude of the commanded relative distance. OC and /? depend on the desired relative motions. In 
the following expression OC and [3 are changing in constant rates 0) a and (Op with the time. 

f a(t) = « 0 + coj 
\P(t) = Pa + a >pt 

The initial values of OC and ft are determined such that the commanded relative distance is equal to the actual 
relative distance at the beginning in order to make the control smooth. 

( 12 ) 


( 11 ) 


r cmdW) = A°r(0) 



Figure 2. Commanded relative trajectory in the target spacecraft orbit frame 


The commanded reference trajectories for several proximity maneuvers are developed. The first proximity 
maneuver is in-track circling, in which the chaser spacecraft circulates the target spacecraft at a circular orbit 
relative to the target with a constant rate. The relative circular orbit is in the target spacecraft’s in-track (x-z) plan. 
The commanded reference trajectory is given by 

( „ , a \ \ 


,(0 = 


' cmd 


'and 


0 0 = 


(0 = 


,„sin(tyt + <9 0 ) 

0 

- r m cos{cot + & 0 ) 
r m cocos(a>t + 0 O ) 
0 


yr m cosin(cot + 0 o ) 


^ — rco 2 sm 1 


(ot + 0 0 ) 


rm~ cos 


{cot + 0 o) 


(13) 


(14) 


(15) 


The magnitude of the commanded relative distance r m is constant for the in-track circling orbit. 0) is the 
constant circling orbit rate. 0 O is determined by Eq.(12). Similarly the cross-track circling is to circulate the target 
spacecraft in the target spacecraft’s cross-track (y-z) plan. The commanded reference trajectory can be expressed as 
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r ow(0 = 


{cot+e,) 


r sm 


y -r m cos{wt + 0 o ) 


( 




0 


\ 


Kmd (0 “ 


r m <ycos(<yr + # 0 ) 

)t 

0 


v r m «sin(<yf + /9 0 ) 


i(<y? + 0 0 ) 


sm( 

v r m or cos (cot + 0 Q ) 


(16) 


(17) 


(18) 


In the next two rendezvous maneuvers the chaser spacecraft approaches the target spacecraft along a spiral 
trajectory in the target spacecraft’s in-track or cross-track plans. For in-track spiral rendezvous the commanded 
reference trajectory is 

f ^(Osin(&>t + <9 0 ) A 

(19) 


( 20 ) 


r cmdi 0 = 

0 



- rJt)cos(cot + 0 () ) y 



' r m (t)sm(cot + 0 O )+ r m (t)cocos(cot + 0 O ) ' 

tmd( t ) = 


0 


v - r m (0 cos(ty t + 0 o )+r m (t)cosm(co t + 0 o \ 


' 2 r m ( f)co cos(ty t + 0 ( 

)-rJl)or sin(ty t + 0 o ) 

Kmd( t ) = 


0 


2r m (t)ojsm(ojt + 0 () 

) + r m ( t)af cos(ty t + 0 O ) 


r nd 

t<t i 


V 

t > t f 


r mi ~ \ r nd ~ r mf N ~ t, )/(tf -tf) t,<t<t 

r.(0 = j 

1 

o > 
1 

1 

t <tj or t > t f 
t,<t<t f 


( 21 ) 


( 22 ) 


(23) 


t i is the starting time of the spiral and t ( is the ending time. r mi is the initial commanded relative distance and 
is the distance at the ending time. For rendezvous r mf is zero. Similarly, the commands for the cross-track 


' mf 

spiral rendezvous are given by 


r ond(t) : 


fcmd(t) : 


r »(Osin(fi>f + 0 o ) 

-r m (t)cos{cot + 8 a )^ 


r m (0 sin(ty t + 0 o )+r m ( t)co cos(ty t + 0 o ) 
- r m (t) cos(ty t + 0 O ) + r m ( t)co sin(ty t + 0 o ) 


(24) 


( 25 ) 
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0 I 

f aml (0 = 2r n Xt)(ocos{at + e 0 )-r n Xt)co 2 sm{cot + e a ) (26) 

2r m (i t)a> sin(<y t + 0 O ) + r m ( t)o 2 COS {cot + do)) 

r m (t) and r m (t) are determined by Eq.(22) and Eq.(23). 

When the control algorithm in Eq.(9) is substituted into Eq.(3), the closed-loop differential equations of the 
proximity relative motion in the inertial frame is given by 


r = K, \ K r ( r cmd - A l‘ r )+ K v (r c md ~ A l^)+ ^ } 


+ M 


R R + r 


■ ■ 3 ■ 


+ 


< c F dc Ai,F dl A'J, 


(27) 




m, m, 


In general, Eq.(27) is a 3 dimensional nonlinear time variant differential equation. In order to measure the 
performance of the algorithm, the difference between the commanded relative trajectory r cmd and the actual relative 
trajectory r is defined as X , which is given by 

(28) 
(29) 


l = r cmd ~Kr 


Eq.(27) can be rewritten in terms of X as 

X + K v X + K r X=f(t) 


with 


/(o =(./;« f y it) m) 

= A', 


, A ‘J, 


f 

R 

X 

R + r 


m t m t 

Lk | 3 

\R + rf j 

m c 

l 


vi 


J 


-(K v+ n ot p ot A?r-2n ol A?r 


The rotation rate vector O) ol of the target spacecraft orbit frame X ol , y ot , Z ot can be written as 


('F = \co xot co yot co zot 




(30) 


(31) 


Matrix Cl ot is given by 


o,= 


0 -CD 


-CD,. 


\ ^ yot 


ZOt 

0 

-CD r . 


yot 


CO 


xot 

0 


(32) 


If the gain matrixes K v and K r are positive definite, and f (t) is zero, it can be proved that X = 0 is a 
globally asymptotically stable equilibrium point of Eq.(29). It means that the actual relative trajectory r will track 
the commanded relative trajectory r cmd and eventually converges to r cmd . Usually f(t) is not zero especially when 

the target spacecraft maneuvers during the chaser spacecraft approaching. However, / (/) can not be unlimited 
because the accelerations due to the target spacecraft thrust and the disturbance forces are limited. If / (/) satisfies 


the bound 


1/(01 ~ /max 

(i = x,y,z) 

(33) 

with 


f =(f f 

J max V a max J 

f f 

y max J zmax / 

(34) 

Then it can be proved that X 

is globally uniformly ultimately bounded as 

(35) 
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as time approaches infinity. The tracking error of the close-loop system with the control algorithm converges and is 
bounded to ||j^ _1 / nax |- The equation shows that increasing the gain matrix K r can reduce the final tracking error. 

Eq.(35) can be used to estimate the final errors between the actual relative trajectory and the commanded relative 
trajectory. The estimate by Eq.(35) is consistent with the simulation results as shown in the next section. 


a°‘t - 

A bt 1 t ~ 


IV. Design Examples 

In order to evaluate the performance of the control algorithm developed in this paper, two rendezvous maneuvers 
are demonstrated. In these two maneuvers, the chaser spacecraft approaches the target spacecraft along a spiral 
trajectory in the target spacecraft’s in-track or cross-track plans. The target spacecraft orbit is 488x528 km with 72 
degrees inclination. The target and chaser spacecraft masses are 600.0 and 400.0 Kg respectively. 

The initial relative distance r( 0) between the target and the chaser spacecraft is set to 10 meters. r( 0) wrt the 
target orbit frame is 

Afr{ 0) = (0 0 -10f (36) 

During the rendezvous the target spacecraft conducts maneuvers by firing the thrusters in all three directions. 
Suppose that the thrust forces in in-track, cross-track and radial direction of the target orbit frame is sine waves and 
given by 

r T x sm(2n:t/t x +0 x )' 

T y sin(2 nt /t y + 0 y ) (37) 

y T z sin(2;zT / t z + 6 Z ) y 

The parameters in the above expression have the following values in the simulations. 

T x = 2.Q(N) T y = 4.0 (N) T =5.0(N) (38) 

t x = 100.0 (j) ? v =60.0(x) t z = 130.0 (j) (39) 

9 X = 20.0 (deg) d v = 80.0 (deg) G z = 120.0 (deg) (40) 

Control gains K r and K v are scalar and given by 

K r = 0.1 (1/s 2 ) A' v =0.1(1/j) (41) 

The parameters in the commanded reference trajectories of the two maneuvers have the following values. 

(0 = 1.0 (deg / s) 9 0 = 0.0 (deg) (42) 

r mi = 10 0 ( m ) r mf = 0-0 O) (43) 

t t = 200.0 (s) t f = 1200.0 ( s ) (44) 

The three components of the chaser commanded thrust forces T c is limited to 8 N. For this two maneuvers, 


./max * s dominated by the target thrust forces and estimated by 

./_-(/ T y Tj/m, (45) 

and 

||^/ 1 / m ax|| = 0.112( W ) (46) 

The actual value of the maximum of |/t| is slightly larger than 0.112 m because f max includes only target thrust 
forces. 

The time histories of the rendezvous maneuver with the chaser spacecraft approaching the target spacecraft along 
a spiral trajectory in the target spacecraft’s in-track plan are shown in Figure 3 to Figure 8. Figure 4 shows that the 

value of the maximum of ||/l||, when it converges, is 0.118 m, which is slightly larger than 0.112 m. The 

commanded thrust forces in in-track direction reaches the limit (8N) during the first 60 seconds. 

Figure 9 to Figure 14 is the time histories of the second rendezvous maneuver. The chaser spacecraft approaches 
the target spacecraft along a spiral trajectory in the target spacecraft’s cross-track plans. The value of the maximum 
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of ||/t|| is 0.118 m also when it converges as shown in Figure 10. Figure 12 and Figure 13 show the chaser 

spacecraft approaches the target spacecraft in the in-track, radial plan and cross-track, radial plan. As we expected, 
the chaser spacecraft approaches the target spacecraft along a spiral trajectory. Figure 14 shows the commanded 
thrust forces calculated by the control algorithm to accomplish the desired rendezvous maneuver. The results show 
the excellent performance of the control system. 


Case 1: In-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 

15 


-10 


-15 I i ' ' i 

0 500 1000 1500 2000 2500 

Time (sec) 

Figure 3. Case 1: Relative Distance 


Case 1: In-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 



Figure 4. Case 1: Relative Distance (zoomed in) 


In-Track 

Cross-Track 

Radial 
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Case 1: In-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 



Case 1: In-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 



Figure 6. Case 1: Relative position in radial and in- 
track plan 
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Relative Velocity (m/s) Relative Dist (m) ^ Radial (m) 


Case 1: In-Track Spiral Case 1: In-Track Spiral 

Spiral Rate = 1 (deg/s) Initial Dist =10 (m) Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 




igure 7. Case 1: Relative position in radial and 
cross-track plan 


Figure 8. Case 1: Thrust forces command 


Case 2: Cross-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 



Case 2: Cross-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 



Figure 10. Case 2: Relative Distance (zoomed in) 


Case 2: Cross-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 



Case 2: Cross-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 



In-Track (m) 


Figure 11. Case 2: Relative Velocity 


Figure 12. Case 2: Relative position in radial and in- 
track plan 
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Case 2: Cross-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 



Cross-Track (m) 


Case 2: Cross-Track Spiral 
Spiral Rate = 1 (deg/s) Initial Dist =10 (m) 



Figure 13. Case 2: Relative position in radial and Figure 14. Case 2: Thrust forces command 

cross-track plan 


V. Conclusion 

A nonlinear trajectory control algorithm of rendezvous with maneuvering target spacecraft is developed in this 
paper. The disturbance forces on the chaser and target spacecraft and the thrust forces on the chaser spacecraft are 
considered in the analysis. The control algorithm developed in this paper uses the relative distance and relative 
velocity between the target and chaser spacecraft as the inputs. A general formula of reference relative trajectory of 
the chaser spacecraft to the target spacecraft is developed and applied to four different proximity maneuvers, which 
are in-track circling, cross-track circling, in-track spiral rendezvous and cross-track spiral rendezvous. The closed- 
loop differential equations of the proximity relative motion with the control algorithm are derived. The tracking 
errors between the commanded relative trajectory and the actual relative trajectory are analyzed. The prediction of 
the tracking errors is obtained. Design examples are provided to show the implementation of the control algorithm. 
The simulation results show that the actual relative trajectory tracks the commanded relative trajectory tightly. The 
predicted tracking errors match those calculated in the simulation results. 
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